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ABSTRACT 

We present new calculations of X-ray polarization from black hole (BH) ac- 
cretion disks in the thermally-dominated state, using a Monte-Carlo ray-tracing 
code in full general relativity. In contrast to many previously published stud- 
ies, our approach allows us to include returning radiation that is deflected by the 
strong-field gravity of the BH and scatters off of the disk before reaching a distant 
observer. Although carrying a relatively small fraction of the total observed flux, 
the scattered radiation tends to be highly polarized and in a direction perpendic- 
ular to the direct radiation. For moderately large spin parameters (a/M > 0.9), 
this scattered returning radiation dominates the polarization signal at energies 
above the thermal peak, giving a net rotation in the polarization angle of 90°. We 
show how these new features of the polarization spectra from BHs in the thermal 
state may be developed into a powerful tool for measuring BH spin and prob- 
ing the gas flow in the innermost disk. In addition to determining the emission 
profile, polarization observations can be used to constrain other properties of the 
system such as BH mass, inclination, and distance. New instruments currently 
under development should be able to exploit this tool in the near future. 

Subject headings: black hole physics - accretion disks - X-rays:binaries 



-2- 



1. INTRODUCTION 

A recent flurry of new mission proposals has renewed interest in X-ray polarization from 
a variety of astrophysical sources. The Gravity and Extreme Magnetism SMEX (GEMS) 



missioro, which has recently been approved for phase A study in the latest round of NASA 
Small Explorer propo sals, should be ab l e to detect a degree of polarization 5 < 1% for a 
flux of a few mCrab (IBlack et al.l 120031 ; iBellazzini et al.l 120061 1 . A similar detector for the 
International X-ray Observatory (IX O) could achieve sensitivity roughly 10 x greater (< 0.1% 
degree of polarization; Jahoda et al. 2007, Costa et al. 2008). A large number of galactic 
and extra-galactic sources are expected to be polarized at the 5 > 1% level, including black 
holes, magnetars, pulsar wind nebulae, and active galactic nuclei. In this paper, we focus 
on accreting stellar-mass black holes (BHs) in the thermal state, which are characterized by 
a broad-band spectrum peaking around 1 keV. The typical level of polarization from these 
sources should be a few percent in the 1 — 10 keV range, depending on BH spin and the 
inclination angle of the accretion disk. 

Symmetry arguments demand that in the flat-space (Newtonian) limit, the observed 
polarization from the disk must be either parallel or perpendicular to the BH/disk rotation 
axis. However, the effects of relativistic beaming, gravitational lensing, and gravito-magnetic 
frame-dragging can combine to give a non-t rivial net rotation to th e integrated polarization 
vect or . Early work exploring these effects (IStark fc Connors! 1 19 771 ; IConnors &: Stark! 119771 : 
Connors et al.lll980l ) showed that they create changes in the angle and degree of polarization 
that are strongest for higher photon energy. The reason for this trend is that the temperature 
in a standard accretion disk increases as the BH is approached, and that is, of course, where 
the relativistic effects are strongest. Since those first efforts, there have been studies of 
the predicted polarization from more general accretion geome tries, including UV and X-ray 



emission from AGN disks (ILaor et al.l Il990l ; iMatt et al.l Il993l ) and "lamp po st" models for 



irrad iating the accreti on disk with a non-t hermal source above the plane (IDovciak et al. 



20041 ). Quite recently, IDovciak et al.l (12008!) investiga ted the effect of atmospheric optical 
depth on the polarization signal, and iLi et al.l (120081 ) applied the original calculations of 
thermal X-ray polarization to the problem of measuring the inclination of the inner accretion 
disk. 

Despite this relatively large body of literature dedicated to X-ray polarization from BH 
accretion disks, to the best of our knowledge, all previous work has modeled the relativistic 
effects by calculating the transfer function along geodesies between the observer and emitter. 
By its very nature, this method precludes the possibility of including the effects of returning 
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radiation — photons emitted from one part of the disk and bent by gravitational len sing so 
that they are absorbed or reflected from another part of the disk f Cunningham|[l~976l ) . Using 
new methodology, in this paper we show the contribution from this returning radiation 
is very important to the observed polarization. The details of our method are described in 



Schnittmanl (120091 ) (hereafter "Paper I"). The most important feature of our approach is that 



the photons are traced from the emitting region in all directions, either returning to the disk, 
scattering through a corona, getting captured by the BH, or eventually rea ching a distant 
observ er. It should be noted that a sort of hybrid technique was used in lAgol fc Krolik 
( 120001 ) to estimate the effects of returning radiation on polarization by situating an observer 
in the disk plane and shooting rays backwards to see how much flux was incident from other 
regions of the disk, but this approach cannot account for multiple scattering events, where 
the returning radiation passes over the BH more than once. 

The inclusion of returning radiation leads to the most important of the new results 
presented in this paper, namely a transition between horizontal- and vertical-oriented polar- 
ization as the photon energy increases. At low energies we reproduce the "Newtonian" result 
of a semi-infinite scattering atmosphere emitting radiation weakly polarize d in a direction 
paral lel to the emission surface, an orientation we call horizontal polarization (I Chandr asekhar 



19601 ). At higher energy, corresponding to the higher temperature of the inner disk, a greater 



fraction of the emitted photons returns to the disk and is then scattered to the observer. 
These scattered photons have a high degree of polarization and are aligned parallel to the 
disk rotation axis (vertical), as projected onto the image plane. At the transition point be- 
tween horizontal and vertical polarization, the relative contributions of direct and reflected 
photons are nearly equal, and no net polarization is observed. 

Since the effects of returning radiation are greatest for photons coming from the inner- 
most regions of the disk, the predicted polarization signature is strongly dependent on the 
behavior of gas near and inside the inner-most stable circular orbit (ISCO). Thus, polar- 
ization observations could be used to measure the spin of the black hole and also constrain 
the dissipation profile of the inner disk. For systems where the BH mass, distance, and disk 
inclination are known, a single polarization observation can constrain the em i ssion profile 



at least well as fitting th e thermal continuum spectrum (jGierririski et al.l l200ll ; iDavis et al. 



20061 ; IShafee et al.ll2006l ). When these system parameters are not known a priori, polariza- 



tion can be used in many cases to measure them simultaneously with the emission profile 
(albeit with lower confidence than for the sources with known priors). 

The structure of this paper is as follows: after briefly describing our methods in Section 
[2j we review in Section [3] some of the classical results of earlier polarization calculations from 
thermal disks, including only the radiation that reaches the observer directly from the disk. 
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In Section HI we include the effects of returning radiation, focusing on new features in the 
polarization spectrum. Sections [5] and [6] contain a detailed discussion of our results and how 
they depend on model parameters, and in Section [7] we conclude. 



METHODOLOGY 



The inner disk of a stellar-mass BH binary c an be well described locally as a plane- 



paral lel, electron scattering-dominated atmosphere (IShakura fc Sunvaev 



1973; Novikov fc Thorne 



19731 ). In the local fluid frame, the classical result of IChandrasekharl (119601 ) applies: the po- 
larization vector of the emitted radiation is oriented parallel to the disk plane and normal to 
the direction of propagation. The degree of polarization varies from zero for photons emitted 
normal to the disk surface up to ~ 12% for an inclination angle of 90°. In addition to the 
polarization effects, the scattering of the outgoing flux causes limb darkening, effectively 
focusing the emitted radiation in the direction normal to the disk surface. Both the degree 
of polarization and the limb dar kening/brighten i ng fa ctors are tabulated as a function of 
emission angle in Table XXIV in IChandrasekharl fll960l ). 



For the local emi ssion spectrum, we assume a diluted black-body s pectrum wit h hard - 



ening factor / = 1.8 (IShimura fc Takaharalll995l ). As demonstrated by iDavis et al.l (120051 ). 
this is a reasonable approximation for disks around stellar-mass black holes because the tem- 
perature is so high that there is little opacity due to anything but Compton scattering and 
free-free absorption. For a locally emitted flux JF, the effective temperature is defined as 



T, 



eS = | — 

a 



1/4 



where a is the Stefan-Boltzmann constant. The spectral intensity of the diluted emission is 
then given by 

I v = j- 4 B u (fT cS ), (2) 

with B v the black-body functi on. For the emitted flux function J-(R), we consider two 
different models: that given by iNovikov fc Thornd (119731 ) for a zero-stress inner boundary 
condition at the ISCO, and a quasi-Newtonian expression ^{R) ~ R~ 3 that extends the 
classical scaling at large radius all the way down to the horizon. Recent global magneto- 
hydrody namic (MHD) sim ulations suggest that reality lies somewhere in between these two 
models ( Noble et ali 2008T). Outside of the ISCO, the gas follows circular, planar geodesic 
orbits as in 



Novikov fc Thornd (119731 ). Inside the ISCO, the fluid continues along plunging 



trajectories with constant energy and angular momentum, resulting in a rapid decline of 
surface density with decreasing radius. 
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As described in Paper I, at each point in the disk, multi-wavelength photon packets are 
emitted in all directions, with appropriate intensity weights and polarizations as a function 
of emission angle. We use a Monte Carlo algorithm with typically 10 5 — 10 6 rays traced 
per radius, with 500 radial bins spaced logarithmically between the horizon and R = 100M. 
Since the photon trajectories originate at the disk, the observer's inclination angle is not 
specified; rather, all photons are followed until they get captured by the BH, return to the 
disk plane, or reach "infinity," in practice a spherical hypersurface at large radius. At that 
point, the photons are binned by observer inclination, effectively simulating all inclinations 
simultaneously. By calculating the angle that each photon makes with the detector sphere, 
we are also able to create an image of the source, much like a classical pin-hole camera. 

Each photon packet contains a complete broad-band spectrum initially described by 
equation (j2J). Since Thomson scattering in the disk atmosphere is frequency-independent 
over the range of relevant X-ray energies (~ 0.1 — 10 keV for thermal emission from a 
solar- mass BH), the photon packet can be described by a single polarization amplitude and 
direction. 



Following IConnors et al.l (119801 ) . we characterize the polarization of each photon packet 
by the normalized Stokes parameters: 

X = Q/I, (3a) 

Y = U/I, (3b) 

and 

5 = (X 2 + Y 2 y'\ (4a) 

^ = itaa-^y/X), (4b) 

where ip is the angle of polarization, ranging from —90° to +90°, and 5 < 1 is the degree of 
polarization, invariant along a geodesic. 

The polarization vector f is a normalized space-like 4-vector perpendicular to the pho- 
ton propagation direction: k • f = and f • f = 1. It is parallel-transported along the 
photon's null geodesic, c onserving the complex-valued Walker-Penrose int egral of moti on k 



( 1 Walker fc Penrose! 1 19701 ). Analogous to the way that Carter's constant (jCarterlll968l ) can 
be used to constrain four- velocity components, k can be used to reconstruct the polarization 
vector at any point along the geodesic. The two or t honor mality conditions stated above, 
along with the real and imaginary parts of k, give a total of four equations for the four 
components of the polarization vector f . When a ray reaches an observer at infinity, the po- 
larization vector is projected onto the detector plane defined by the basis vectors ei and e 2 to 
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give the polarization angle: f = [0, cos^ ei + sinip (in our convention, e2 is parallel to the 
projected rotation axis of the BH and disk). From 5 and i/j, as measured by the detector, we 
can reconstruct the normalized Stokes parameters X and Y via equation (J3J). The observed 
spectrum I v is propagated along the geodesic path as well, containing within it the emitted 
spectrum, modified by all relativistic effects as well as any energy change that may occur 
when returning radiation scatters off disk electrons. I u provides the absolute normalization 
that, when applied to X and Y allows Q u and U u to be recovered for each photon bundle 
reaching the detector. Integrating over all photon bundles gives the polarization degree 8 V 
and angle ip v as a function of energy (see Paper I for details). These are ultimately the 
quantities that we wish to be able to measure and use to probe the accretion geometry of 
the BH. 



3. DIRECT RADIATION 

In Figure [TJ we show a simulated image of a Novikov-Thorne accretion disk around a 
black hole with spin parameter a/M = 0.9 and luminosity O.ILe, corresponding to a disk 
whose X-ray spectrum peaks around 1 keV for a BH mass of 1OM . With the observer at 
an inclination of i — 75°, significant relativistic effects are clearly apparent. The increased 
intensity on the left side of the disk is due to special relativistic beaming of the gas moving 
towards the observer, and the general relativistic light bending makes the far side of the disk 
appear warped and bent up towards the observer. Superposed on top of the intensity map is 
the polarization signature, represented by small black vectors whose lengths are proportional 
to the degree of polarization observed from that local patch of t he disk. Far from the black 



hole, the polarization is essentially given by the classical result of lChandrasekharl (119601 ) for a 
scattering-dominated atmosphere: horizontal orientation with 5 ~ 4% when % = 75°; nearer 
to the black hole, a variety of relativistic effects alter the polarization. 



While the 



undamental calcula tion of polarization from a thermal BH accretion disk is 



not at all new (I Connors et al.l Il980l ) , to the best of our knowledge this is the first published 
polarization map as seen projected in the image plane. The usefulness of such a map is 
primarily as a tool to understand how relativistic effects determine the integrated polarization 
at different energies. 

The two most prominent relativistic effects are gravitational lensing and special relativis- 
tic beaming, both lowering the net level of polarization seen by the observer. Gravitational 
lensing causes the far side of the disk to appear warped up towards the observer, and thus 
have a smaller effective inclination. Relativistic beaming causes photons emitted normal to 
the disk plane in the fluid frame to travel forward in the direction of the local orbital motion 
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when seen by a distant observer; the result is a smaller effective inclination and thus degree of 
polarization. On the other hand, photons emitted at high inclination angle in the local fluid 
frame but against the direction of orbital motion are ultimately seen by observers located 
at lower inclination, and therefore raise the polarization in that direction. Naturally, these 
relativistic effects are most important close to the black hole, where the gas is also hottest 
and the photons have the highest energies. All these effects are clearly visible in Figure [TJ 
which shows a smaller degree of polarization where the beaming is greatest (yellow region 
of high intensity in the left of the image) and the lensing is strongest (just above the center 
of the image). At the same time, the gas moving away from the observer on the right of the 
image has an enhanced level of polarization because the observer sees photons emitted at a 
larger inclination in the fluid frame. 

Integrating over the entire disk, we can calculate the angle and degree of polariza- 
ti on as a function of en ergy. Plotted in Figure El these results essentially reproduce those 



of I Connors et al.l (119801 ). who performed a similar calculation using semi-analytic transfer 
functions from the observer to the disk. Again we see that at low energies (dominated by 
emission from large disk radii), the polarization is given by the Chandrasekhar result, and is 
a function only of the disk inclination angle. At higher energies, we begin to probe the inner 
regions of the accretion disk and see the relativistic effects described above that lower the 
total degree of polarization. In addition to reducing the polarization amplitude, the strong 
gravity near the black hole also rotates the angle of polarization due to the parallel transport: 
since the polarization vector f must remain perpendicular to the photon momentum k, as the 
geodesic path bends around the BH, the polarization angle must also rotate. For example, 
the high-energy photons beamed towards the observer in Figured] are also bent by the BH's 
gravity, initially moving to the left before curving back to the observer, thus producing the 
rotated polarization vectors seen just to the left of the image center (the same lensing effect 
rotates the red-shifted photons from the right side of the disk, but in an opposite direction). 

The angle of polarization as a function of energy is plotted in the lower panel of Figure [2j 
For an inclined accretion disk with the black hole axis projected onto the vertical, or y-axis, 
ip = corresponds to horizontal polarization parallel to the disk surface. From equation (jlj), 
we see there is a point symmetry in the definition of ip, giving if> = if> ± 180°. As can be 
seen in the simulated image of Figure [TJ the individual polarization vectors are rotated in the 
inner disk, due largely to gravitational lensing and, to a lesser degree, frame-dragging around 
a spinning black hole. While individual photons can experience significant rotation, the net 
result is a modest rotation of the total observed polarization in the clockwise direction, 
giving if) < at higher energies. From Figure [21 it is clear this effect is more pronounced for 
small inclinations, where the low amplitude of polarization makes it easier to "overcome" 
the classical result with additional relativistic effects. Similarly, the polarization rotation 
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is greater for more rapidly spinning black holes, where the accretion disk extends closer in 
towards the horizon and probes a stronger gravitational field. In this way, it has previously 
been proposed that the polarization degree and angle as a function of energy could be used 



to infer the spin of t he black hole in the thermal state (IConnors et al.lll980l ; lLaor et al.lll990 



Dovciak et all 120081 ). 



4. RETURNING RADIATION: QUALITATIVE DISCUSSION 

When returning radiation is included, although little changes in terms of the total 
observed spectrum, the polarization picture (Fig. [3]) changes significantly — in much of the 
inner disk, the integrated polarization rotates by 90°, even though none of the model's 
physical parameters has been changed at all! 



This effect can be understood qualitatively in very simple fashion (see also lAgol fc Krolik 



(120001 ) ) . For most reasonable stellar- ma ss BH accretion disk models, the opacity in the inne r 



disk is dominated by electron scattering (IShakura fc Sunyaevlll973l ; iNovikov fc Thorndl 19731 ) 



so returning radiation in the ~ 1 — 10 keV energy range should scatter off the disk with 
negligible absorption. Detailed atmosphere calculations (S. Davis, private communication) 
show that the photospheric region is strongly scattering dominated for R < 50M whenever 
-^/-^Edd ^ 0.01. If these calculations were altered to take into account magnetic contributions 



to ve rtical support, which should be substantial at these altitudes in this regime (IHirose et al. 



20091 ). absoroption opacity would likely be even weaker (however, when the central mass is 
much larger, as in an AGN, the inner disk temperature is much lower, giving significantly 
greater absorptive opacity in the X-ray band). Moreover, the typical photon scatters only 
once or twice before permanently departing the disk. Its outgoing polarization is therefore 
very sensitive to the electron scattering cross section's polarization-dependence: 



dtt) 



^■f/l 2 , (5) 



pol 



where r is the classical electron radius. In other words, scattering is strong only when the 
polarization direction changes little. This is possible for both senses of linear polarization 
when the change in photon direction is small, but when the scattering is nearly perpendic- 
ular, the new photon wave-vector becomes nearly parallel to one of the initial polarization 
directions. Thus, for nearly-perpendicular scattering, the outgoing light can be close to 100% 
polarized. 

For observers at high inclination angles, such as in Figure El returning radiation photons 
initially emitted from the far side of the disk (top of the image) are reflected off the near 
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(bottom) side with a relatively small scattering angle, maintaining a moderate horizontal 
polarization as in Figure [TJ On the other hand, photons emitted from the left side of the disk 
can be bent back to the right side (or vice versa), and then scatter at roughly 90° to reach 
the observer, thereby aquiring a large vertical polarization component. Although relatively 
small in total flux, this latter contribution can dominate the polarization because it is so 
strongly polarized. A flux component of only 10%, when 100% polarized, contributes more 
to the polarized flux than does a 90% contribution that is only 5% polarized. 

Because high-energy photons from the hotter inner parts of the disk experience stronger 
gravitational deflection, they are more likely to return to the disk than the low-energy pho- 
tons emitted at larger radii. We find that, for each value of the spin parameter, there is 
a characteristic "transition radius," within which the returning radiation dominates and 
produces net vertical polarization. Outside this point, the direct radiation dominates and 
produces horizontal polarization. This transition radius is in fact only weakly dependent on 
a/M, ranging from R trans « 7M for a/M = to R tvSLlls « 5M for a/M = 0.998. The location 
and shape of the observed polarization swing can be used to infer the radial temperature 
profile near the transition radius (see Sec. [6] below for more details). While the majority 
of the returning obviously originates from the inner-most disk, we find that most of it also 
scatters inside R < 10M as well (where the atmosphere opacity is completely dominated by 
electron scattering), due to the additional focusing effects of gravitational lensing. 



5. RETURNING RADIATION: QUANTITATIVE RESULTS 

Polarization maps like that shown in Figure [3] provide a useful picture for understanding 
the physical origin of the local polarization features, but are less useful in providing quantita- 
tive predictions for future observations. Since vertically- and horizontally-polarized photons 
add up to give zero net polarization, it is not clear from Figure [3J alone what the integrated 
polarization signal should be. Furthermore, Figure [3J shows only the energy-integrated polar- 
ization at each point in the disk, masking the important energy dependence of polarization 
seen in Figure [2j 

In Figures HH6] we show the spectral intensity, degree, and angle of polarization as a 
function of photon energy for a range of BH spins and inclinations. In each frame, we plot 
the total flux (solid curves), as well as the direct flux from the disk (dotted curves) and 
the returning radiation (dashed curves). Beginning with the spectral intensity, we see the 
standard broad thermal peak characteristic of BHs in the "high soft" state. For BH masses 
around 1OM and luminosities of ~ O.lZvEdd, the thermal spectrum peaks around 1 keV. For 
a Novikov-Thorne (NT) emission profile, the direct radiation from a non-spinning BH (peak 
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emissivity at radius R ~ 9M) dominates over the returning radiation by a factor of roughly 
100. As the spin increases and the peak emission region moves closer to the horizon, the 
relative fraction of observed flux that comes from returning radiation increases to ~ 5% for 
a/M = 0.9 and ~ 20% for a/M = 0.998. Naturally, the relative contribution from returning 
radiation increases with energy because the highest energy photons come from the innermost 
regions of the disk. In fact, for a/M = 0.998, the returning radiation actually dominates the 
spectrum above ~ 10 keV. 

For all spins, the fraction of returning radiation increases with inclination for two rea- 
sons: forward/backward scattering has a larger differential cross-section than right-angle 
scattering, so returning photons, which characteristically have large angles of incidence, are 
reflected at large angles from the axis and preferentially reach observers with high incli- 
nation; secondly, the direct radiation suffers limb-darkening effects at high inclination due 
both to the scattering atmosphere and also the diminished solid angle of the disk, while both 
effects affect the reflected radiation to a lesser degree (see Paper I for a detailed calculation). 
For the examples shown in Figures HHfH the fraction of observed flux coming from returning 
radiation is larger for i = 75° than i = 45° by a factor of about two. 

At low energies (below the thermal peak), the flux comes predominantly from the 
outer regions o f the d isk, where relativistic effects are negligible, so the classic results of 
Chandrasekhar Jl960h are reproduced horizontal polarization with amplitude of a few per- 



cent. The degree of polarization at low energies increases with observer inclination, ranging 
from 5 = 1% for i = 45° to 5 = 3.5% for i = 75°. Although the free-free absorption opacity 
in the disk increases below ~ 1 keV, which should reduce the net polarization, we have 
ignored this effect for two reasons. First, large uncertainties in the vertical structure of the 
disk make it difficult to arrive at a quantitative prediction for the absorption profile, which 
is strongly dependent on local temperature and density. Second, most proposed X-ray po- 



larization missions will be sensitive only to photons above 1 keV (jBellazzini et al.ll2006l ). so 
we choose to focus on this energy range. 

Above the thermal peak, we see the relativistic effects described in the previous Sections 
begin to dominate the polarization signal. The direct radiation (dotted curves) undergoes 
lensing and beaming, thereby reducing the polarization amplitude and rotating the angle. 
The polarization properties of the returning radiation (dashed curves) are roughly constant 
in the 1 — 10 keV range, but as the relative flux from returning radiation increases, we see 
the total polarization makes a transition from following the direct behavior to following the 
returning behavior. In the process of making this transition from horizontal polarization at 
low energies to vertical polarization at high energies, the degree of polarization goes through 
a minimum as the two contributions cancel each other. This transition is sharpest for lower 
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inclinations, where there is a smaller degree of horizontal polarization to overcome. Yet 
even for large inclination angles, when the BH spin is high enough, we still see a dramatic 
transition from horizontal to vertical polarization close to the thermal peak, caused by the 
enhanced level of returning radiation from the innermost disk. 

Due to this dependence on the relative contribution from returning radiation, if the 
emissivity profile is NT, the energy at which the polarization transition takes place may be 
a powerful diagnostic of BH spin. In Figure [7] we plot the polarization degree and angle as 
a function of energy for a range of spin parameters. In all cases, the BH mass is 1OM , 
i = 75°, and the luminosity is O.lLEdd (due to differences in accretion efficiency, we use a 
fixed luminosity, not mass accretion rate, when comparing different spins). For inclinations of 
i = 45° or 60°, this same dependence on spin can be seen by comparing the top and middle 
rows of Figures 0H2J With even moderate energy resolution and polarization sensitivity 
of < 1%, observations of BHs in the thermal state should be able to distinguish between 
Schwarzschild and extremal Kerr with high significance — if one can safely assume the NT 
surface brightness profile. 

However, recent relativistic MHD simulations suggest that, contrary to the NT assump- 
tion of zero stress and zero dissipation inside the ISCO, there are in fact significant stress 
and dissipation in the plunging region, leading to steadil y increasing; surfac e brightness (as 



measured in the local fluid frame) with decreasing radius fjNoble et al.ll2008l ). For a simplis- 
tic treatment of this flux, we apply a quasi-Newtonian emissivity profile with ^{R) ~ i?~ 3 
all the way to the horizon, regardless of the BH spin or ISCO location. The local spectral 
distribution is still given by equations (|T]{2j) . The resulting polarization signals are plotted 
in Figure [BJ showing a clear reduction in the sensitivity to spin parameter. While there is 
still some variation between the different curves, it would likely be impossible to measure 
a/M with a first-generation X-ray polarization mission (see below, Sec. E]). 

The transition energy in the polarization angle is a strong function of the luminosity, 
which in turn is a function of the temperature of the inner disk. Since the initial polarization 
and subsequent scattering-induced polarization are independent of frequency at these pho- 
ton energies, we can think of the luminosity dependence simply as a rescaling of the entire 
polarization spectrum with inner disk temperature. This scaling is illustrated in Figure 
which plots the polarization signature for a range of luminosities, in all cases using a NT flux 
profile with a diluted thermal spectrum and BH mass 1OM , for a/M = (solid curves) and 
a/M = 0.9 (dashed curves). Since most galactic BHs are observed over a wide range in lu- 
minosities, even within the thermal state, multiple polarization observations could give even 
better constraints on the value of the spin parameter. By reproducing the same polarization 
at a shifted energy, such a measurement would also help confirm that the features in the 
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polarization spectrum are in fact coming from relativistic effects in the inner disk. Again, 
if the underlying emission profile increases continuously down to the horizon (see Fig. [TO]) , 
we lose some sensitivity in determining the spin, but gain sensitivity to L/L-^^ because the 
dependence on that variable sharpens. 



6. PARAMETER ESTIMATION 

As shown in the previous section, spectropolarimetry measurements of BHs in the ther- 
mal state will give a powerful new way of measuring BH spin, if the underlying emission 
follows a Novikov-Thorne profile with zero dissipation inside of the ISCO. More generally, 
polarization observations will probe the emissivity profile directly, offering a new way to 
measure the temperature profile as a function of the geometrized radius (R/M). Well out- 
side the ISCO, conservation laws constrain all models of time-steady disks to follow the NT 
form; however, nearer the ISCO and throughout the plunging region inside it, m ore variation 
i s pos sible. Motivated by global MHD simulations such as those presented in iNoble et al. 



(120081 ) and iBeckwith et al.l (120081 ). we describe the range of possible dissipation profiles by 
power-laws in radius in the plunging region matched smoothly to the NT profile just outside 
the ISCO. Figure [11] shows a sample of emission profiles for a/M = 0.9 and a range of 
power-law index a. In the limit a — * — oo, we reproduce NT; a = 3 is the pseudo-Keplerian 
case shown above in Figures [H] and [TUJ Also shown for comparison is the dissipation profile 
calculated by HARM3D for th e case of a/M = .9 and a disk aspect ratio H/R ~ 0.1, which 



corresponds to a ~ 0.5 — 1 (INoble et al.l 120081 ) . 



To estimate the ability of polarization observations to accurately determine a/M and 
a for a typical galactic BH binary, we consider two theoretical X-ray polarimeters: a "first- 
generation" instrument with modest energy range 1 — 10 keV, resolution AE/E « 1, and 
minimum polarization sensitivity of 5 > 1% with la confidence. A "next-generation" detec- 
tor is characterized by a broader energy range 0.1 — 100 keV, better resolution AE/E « 0.1, 
and ten times the collecting area, giving 5 > 0.3% for the same source. We then select a few 
sample "target" models with different emission profiles and generate simulated polarization 
data for these cases using the ray-tracing code described above. Scanning over a range of 
model parameters, we calculate a large number of "template" polarization spectra to be 
matched against the targets. 

For each point in parameter space, we calculate the overlap between a template spec- 
trum and the target spectrum by calculating y 2 for the specified detector sensitivity and 
source luminosity. Rather than attempting to combine the absolute luminosity, distance, 
detector response function and collecting area into an error estimate for a generic spectro- 
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polarimeter, we simply describe the instrument and source together in terms of the minimum 
error achievable in a measurement of 5 for a single energy bin. For a broad-band detector 
with constant sensitivity at all energies, the minimum error will correspond to the peak of 
the observed intensity spectrum, i.e. the bin with the most photon counts. The error in 
any other bin simply scales like the square root of the counts in that bin. We use for our 
observables the Stokes parameters Q u and U u , both normalized to the total observed flux 
hot — f Ivdv. For a peak sensitivity (minimum error) of £ min , the error on each value of Q u 
and U v is given by 

AQ u = AU u = S min I v ^^. (6) 

With the target values of Q u and U v and their associated measurement errors at each energy, 
it is straightforward to calculate the x 2 fit for the match between the target "observation" 
and any given template spectrum. By generating multiple templates to cover the parameter 
space, we can estimate the significance with which the model parameters can be determined 
from future data. 

In practice, this search through parameter space would be prohibitively expensive com- 
putationally if each template were constructed using the full Monte Carlo calculation (typi- 
cally ~ 10 8 rays traced per run). Fortunately, we are able to utilize transfer functions similar 
to those used in c alculating the relativistic effects on broad-band continuum spectra (e.g. 



Davis et al.l (120051 )). The full details are described in Paper I, but a short summary of the 
method is as follows: for each value of the spin parameter (roughly 20 points spaced evenly 
in log(l — a/M)), a single complete Monte Carlo simulation is carried out as above, with 
the emission spectrum at each radius given by a delta-function in energy. The initial polar- 
ization is still given by a scattering-dominated atmosphere, and the returning radiation is 
treated as before, typically adding a strong vertical component to the observed polarization. 
When each ray reaches an observer at infinity, the Stokes parameters I u , Q u , and U v are 
recorded, each looking essentially like a delta function shifted in energy and amplitude by 
the relativistic effects of the ray-tracing. 

Sorting all the rays by observer inclination angle, we construct a transfer function 
g(a/M) that takes the input emission profile I„ m (R) and generates the observed polarization 
in terms of I° bs (i), Q° bs {i), and U° hs (i). For a given spin, a single transfer function can be 
used to generate polarization spectra for any set of model parameters M, L/L Edd , and a, 
all of which combine to give I v ^ l {R). The transfer function can also be used to calculate 
the transition radius -Rtrans described above in Section HI For emission profiles with steep 
gradients around -Rtrans (as would be produced by large values of a/M or a; see Figs. |8]and 
[TO]) , the observed polarization spectrum displays a sharp swing around the photon energy 
corresponding to the disk temperature at that radius: -Et rans ~ 3T(i? trans ). On the other 
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hand, moderate spin systems with NT emission profiles have a shallow temperature gradient 
around -Rtrans; giving a broad swing in the polarization angle (Figs. [7] and [9]), or none at all 
in low-spin cases where there is simply not enough emission from inside -Rtrans to overcome 
the horizontal contribution from emission outside of -Rtrans- 

For reasonably high resolution spectra, we find the transfer function method gives a 
speed-up factor of better than 10,000 compared to the direct Monte Carlo calculation. This 
method is particularly straightforward in the thermal state, where the scattering cross sec- 
tions and polarization transport are essentially independent of energy. However, even when 
including inverse-Compton coronal scattering, the transfer function method is still applica- 
ble, with the added dimension that the transfer function is dependent on the seed photon 
energy and accretion geometry (again, see Paper I for details). 

In Figure d21 we show quality-of-fit contours for matching the target parameters of three 
different models: two NT radial profiles with a/M = and a/M = 0.998, and one case with 
additional emissivity in the plunging region with a = 1 and a/M = 0.97. All three cases 
have M = 10M Q , L/L Edd = 0.1, and i = 75°. The top row of figures corresponds to cases 
where the disk inclination and Eddington-scaled luminosity (i.e. the BH mass and distance) 
are known a priori, and the target data is fit by scanning over a and a. The middle row 
repeats this calculation, but assumes prior knowledge only of the disk inclination, and then 
for each value of a and a minimizes \ 2 over L/L-^dd- Finally, the bottom row assumes no 
prior knowledge about the system and minimizes x 2 with respect to L/L Edd and inclination 
i. In each frame, the target parameters are marked with an 'X' and the contours show 
confidence intervals in significance units: < lcr (white), 2a (blue), 3a (purple), 4cr (red), 5a 
(orange), and > 5a (yellow). Here a is given by a standard \ 2 distribution: a = V2u, with 
v the number of data points minus the number of free parameters. For large values of 
these confidence intervals correspond very closely to those of a normal distribution: (68%, 
95.5%, 99.7%, etc.). 

As seen in Figure [9], the Schwarzschild case closely resembles the Chandrasekhar limit 
of horizontal polarization with amplitude of a few percent. When the emission cuts off at the 
ISCO, relativistic effects play a much smaller role in rotating the polarization of the direct 
radiation. Furthermore, the fraction of returning radiation is much smaller, minimizing the 
contribution from vertically polarized photons. The only way to generate such a polariza- 
tion spectum is by producing the observed flux at large radii, from where the photons can 
propagate through nearly flat space to the distant observer, arriving with their polariza- 
tion essentially unchanged. Thus if the observed polarization is relatively strong (> 1%), 
horizontally oriented, and constant with energy, then the emission likely comes from large 
radii, in turn constraining the spin to be relatively small and strongly limiting the amount 
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of dissipation inside the ISCO. This is evident from the upper-left frame in Figure [T2l where 
a polarization observation of a Schwarzschild BH with a NT dissipation profile rules out the 
possibility of a high value for a/M or a. 

On the other hand, when matching the polarization from a near-maximal spin BH with 
a/M = 0.998 (upper-right frame of Fig. [T2|) . any system with a NT emission profile and spin 
a/M < 0.97 is ruled out at the 3a confidence level. Yet as described above in Section [5], a 
non-spinning BH can appear to be rapidly spinning if there is significant emission inside the 
ISCO. Additionally, since the plunging region is so small for large values of a/M, templates 
with a range of a (i.e. the emission profile inside the ISCO) are nearly degenerate, as seen 
in Figure [8j These systems are characterized by a sharp transition in the polarization: it 
changes from strong, horizontal polarization at low energies to strong, vertical polarization 
above the thermal peak. 

For intermediate cases with moderately high spin and a [(0.9 < a/M < 0.99); (—2 < 
a < 2)], there is a broad minimum in the polarization amplitude between 1 and 10 keV 
(see, e.g., the red and orange curves in Fig. [7]). If the inclination is known to be relatively 
high, as in our test cases, then this broad minimum can still constrain a/M and a (or at 
least some combination of the two parameters) reasonably well, as seen in the upper-middle 
plot of Figure [12j However, if we have no prior knowledge of the disk inclination, we could 
reproduce such a feature by assuming a low-inclination system, whose light would have little 
or no polarization near the thermal peak regardless of the emission profile. Therefore, as 
shown in the lower-middle panel of Figure [T2"} when minimizing x 2 over all model parameters, 
it is always possible to get a decent fit to an unpolarized source, and the emission profile is 
poorly constrained. At the same time, if a non-zero polarization is detected, it can give a 
measure not only of the emission profile, but also of the disk inclination. 

There are also a number of BH binary systems for which the inclination is known 
relatively well, but the BH mass M and/or distance D are not known. This is the situation 
illustrated in the middle row of Figure [12j As those panels show, substantial constraining 
power is retained even without knowledge of L/Lgdd (and as the bottom panels show, even 
when the inclination angle is also unknown, spectropolarimetry still provides some limit on 
the possible range of a and a/M). The reason is that with a single spectropolarimetry 
observation, we can measure, using the shape and energy of the polarization transition, the 
temperature profile T(r) in the inner disk as a function of radius in geometric units r = R/M. 
From this temperature profile, we can constrain the system's ratio of bolometric luminosity 
to mass through the relation: 



Using the net flux gives Lj D 2 , so that we have two constraint equations for three unknowns. 




(7) 
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If either the BH mass or the distance can be measured by some other means, all three 
variables (L, M, and D) will be known. 

Viewed in this way, spectropolarimetry measurements can be seen to have significant 
advantages relative to pure spectral measurements. When the inclination, mass, and distance 
are all known, detailed continuum-fitting can be used to measure the temperature profile in 
the inner disk, giving constraints on a and a/M that are comparable to what is shown in the 
top row of Figure [T2J However, when one or more of these quantities (inclination, distance, 
or mass) is not known, the continuum-fitting method becomes less able to constrain the 
radial emission profile and thus the BH spin. 

At the qualitative level, one can understand the greater power of spectropolarimetry as 
simply the result of having more observables: 5(E ohs ) and ip(E ohs ) in addition to F(E ohs ), 
giving more constraints on the model parameters. However, there are more specific reasons as 
well. On the one hand, inferring parameters from continuum spectra can be very ambiguous 
because there are multiple ways to reproduce an observed spectrum by varying M, D, i, and 
L/L-Edd- On the other hand, especially for high-inclination systems, the polarization tends 
to constrain individual parameters more directly: for example, its magnitude at low energies 
(where relativistic effects are weak) directly relates to the inclination. At higher energies, the 
distinctive swing in polarization angle is a direct probe of extreme gravitational lensing and 
returning radiation, giving a sensitive indicator of strongly relativistic effects in the inner 
disk. 

Not surprisingly, a next-generation polarimeter with roughly a factor of ten improvement 
in energy bandwidth, resolution, and collecting area can do a much better job at measuring 
all these system parameters. Figure [13] shows the same confidence contours as Figure [121 
now for the next-generation detector. Almost all degeneracy has been removed, and even 
when we have no prior knowledge of the binary parameters, a/M, a, i, L/M 2 , and M/D 
can be determined with high precision. For the three cases considered in the bottom row of 
Figure [T3| the inclination is recovered within ~ 5° and the accretion rate within ~ 20%: 



In fact, at these levels of polarization sensitivity, the ability to measure the intrinsic 
properties of the system likely becomes more dependent on the accuracy of the underlying 
emission model, rather than on the signal-to-noise of the observation. For example, the 
simplified form of the emission spectrum may be modified by a more careful treatment of 
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° < z < 80° 




and 
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the disk atmosphere (IDavis et al.ll2005l ). Additionall y, the low-energy polarization may be 
strongly modified or reduced by absorption in the disk (ILaor et al.lll990l ). Faraday rotation as 
the photons pass through turbulent magnetic fields (S. Davis, private communication), and 
even a small amount of inverse- Compton scattering that is sometimes seen in the thermal 
state. 



7. DISCUSSION 

We have presented here the first results of a new Monte Carlo ray-tracing code for 
calculating the X-ray polarization from black holes. This code, described in detail in Paper 
I, is most notable for its emitter-to-observer paradigm of radiation transport, which allows for 
the inclusion of returning radiation and electron scattering in a completely general accretion 
geometry. In this paper we focus on polarization signatures of the thermal state in stellar- 
mass black holes, a condition in which the disk is optically thick and geometrically thin, and 
its opacity is scattering-dominated. The emitted radiation has a diluted thermal spectrum 
and is weakly polarized parallel to the disk surface. The integrated polarization spectrum 
seen by a distant observer contains distinct energy-dependent features, as the high-energy 
photons from the inner disk are modified by relativistic effects such as Doppler boosting and 
gravitational lensing near the BH. 

For radiation originating very close to the BH, the most important relativistic effect is 
the strong gravitational lensing that causes the photons to get bent back onto the disk and 
scatter towards the observer. This scattering can induce very high levels of polarization, 
especially at large inclination angles, and leads to a distinct transition in the polarization 
angle from horizontal at low energies to vertical above the thermal peak. Observing such a 
swing in the polarization angle would give the most direct evidence to date for the extreme 
relativistic light bending predicted around black holes. Furthermore, by measuring the loca- 
tion and shape of this transition, we will be able to constrain the temperature profile of the 
inner disk. If we assume a NT disk with zero emission from inside the ISCO, the polarization 
transition energy gives a direct measurement of the BH spin (see Fig. [7]). Alternatively, if 
we relax this assumption, the polarization can constrain models for the dissipation profile 
and provide compelling evidence for strong-field gravitational effects. 

For BH systems where the mass, inclination, and distance are known from other obser- 
vations, polarization measurements would be comparable in power t o the continuum fitting 



method (jGierliriski et all 1200 ll ; IDavis et all 120061 ; IShafee et al.1 120061 ) for the purpose of in- 



ferring the emissivity profile of the inner disk. However, when we lack knowledge of any one 
of these priors, polarization provides a significantly stronger tool than continuum fitting for 
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constraining both the shape of the emissivity profile and the unknown parameter (s). 

In this paper, we have presented the first results from our new black hole X-ray polar- 
ization analysis code, corresponding to thermal emission from a geometrically thin accretion 
disk. In future work, we will extend its application to cases in which a hot corona par- 
tially covers a cooler disk (as may be the case in AGN) and to other accretion geometries, 
appropriate to other spectral states of Galactic black hole binaries. 
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Fig. 1. — Ray-traced image of direct radiation from a thermal disk. The observer is located 
at an inclination of 75° relative to the BH and disk rotation axis, with the gas on the left 
side of the disk moving towards the observer, which causes the characteristic increase in 
intensity due to relativistic beaming. The black hole has spin a/M = 0.9, mass M = 1OM0, 
and is accreting at 10% of the Eddington limit with a Novikov-Thorne zero-stress emissivity 
profile, giving peak temperatures around 1 keV. The observed intensity is color-coded on a 
logarithmic scale and the energy-integrated polarization vectors are projected onto the image 
plane with lengths proportional to the degree of polarization. 
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Fig. 2. — Polarization degree and angle from a thermal disk as a function of observed photon 
energy, including only direct radiation. The BH parameters are as in Fig. [1] (dashed curves) 
and also include identical calculations corresponding to a non-spinning BH (solid curves). 
The angle of polarization is measured with respect to the horizontal axis in the image plane 
with point symmetry through the origin: if; — ip — 180°. 
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Fig. 3. — Ray-traced image of radiation from a thermal disk, as in Fig. [TJ but here including 
returning radiation. Photons emitted from the inner disk get bent by the BH and scatter off 
the opposite side of the disk towards the distant observer. 
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Fig. 4. — Intensity spectrum and polarization degree and angle from a thermal disk, including 
direct and returning radiation, for a BH with a/M = 0, M = 1OM , L/L^dd = 0.1, and a 
Novikov-Thorne emission profile. In the left column we plot the observed flux, in the middle 
the degree of polarization, and on the right the angle of polarization. In each plot, the flux 
is divided into contributions from the direct (dotted curves), reflected (dashed curves), and 
total (solid curves). From top to bottom, the observer inclination is 45°, 60°, and 75°. 
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Fig. 5. — Intensity spectrum and polarization degree and angle from a thermal disk, as in 
Figure H, but for a Kerr BH with a/M = 0.9. 
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Fig. 6. — Intensity spectrum and polarization degree and angle from a thermal disk, as in 
Figure H but for a Kerr BH with a/M = 0.998. 




Fig. 7. — Polarization degree and angle for a range of BH spin parameters. All systems have 
inclination i = 75°, BH mass 1OM , luminosity L/L-^dd = 0.1, and Novikov-Thorne radial 
emission profiles. 
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Fig. 8. — Same as Figure [7] but for a power-law emissivity profile ^(R) 
to the horizon. 
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Fig. 9. — Polarization degree and angle for a range of luminosities for a/M = (solid curves) 
and a/M = 0.9 (dashed curves). All systems have inclination i = 75°, BH mass 1OM , and 
Novikov-Thorne radial emission profiles. 
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Fig. 10. — Same as Figure M but for a power-law emissivity profile ^{R) ~ R 3 all the way 
to the horizon. 



100.0 



CD 

CD 
1— 

O) 

CD 

T3 

o 
cc 

N 
'i_ 
CT5 

O 
Q. 



10.0 



1.0 



0.1 



0.1 



Edd 0.03 




1.0 

E o b s (^ V ) 



10.0 



CD 
73 

_CD 

cd 

o 

'"5 

N 

o 

Q. 



-100 



0.1 



Edd 0.03 




1.0 

E o b s ( keV ) 



10.0 



-31 - 



Fig. 11. — Local flux (as measured in the fluid frame) versus radius for a/M = 0.9. In the 
plunging region, the flux is parameterized by a power law with index a, normalized to match 
the NT profile (solid curve) just outside the ISCO. Also shown is the emissivit y calculated 
by HA RM3D (dotted curve), a 3-dimensional global MHD code in full relativity (INoble et al. 
2008h . 
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Fig. 12. — Contour plots of fit quality in matching emissivity profiles to simulated data for 
a first-generation X-ray polarimeter. The left column corresponds to simulated data from a 
non-spinning Schwarzschild BH with a NT emissivity profile (a = — oo), the center column 
has a/M = 0.97 and a — 1, and the right column has a/M = 0.998 and a = — oo. The top 
row shows the fit quality regions assuming the luminosity in Eddington units and the disk 
inclination are both known; the middle row assumes we know the disk inclination but not 
the accretion rate, and the bottom row assumes we know neither. The color coding is < la 
(white), 2cr (blue), 3a (purple), 4cr (red), 5a (orange), and > 5a (yellow), where a is the 
variance of the x 2 ^) distribution, as defined in the text. 
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Fig. 13. — Contour plots of fit quality in matching emissivity profiles to simulated data for 
a next-generation X-ray polarimeter with broader energy range and larger collecting area. 
The different frames and their color codes are the same as in Figure [T2l 
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